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We investigate the diffusive motion of an overdamped classical particle in a f D random poten- 
tial using the mean first-passage time formalism and demonstrate the efficiency of this method in 
the investigation of the large-time dynamics of the particle. We determine the £o<j-tfme diffusion 
{{x 2 (t)) th } dis = Aln 13 (t/t r ) and relate the prefactor A, the relaxation time t r , and the exponent 
(3 to the details of the (generally non-gaussian) long-range correlated potential. Calculating the 
moments {{O th ) dis °f the first-passage time distribution P(t), we reconstruct the large time dis- 
tribution function itself and draw attention to the phenomenon of intermittency. The results can 
be easily interpreted in terms of the decay of metastable trapped states. In addition, we present a 
simple derivation of the mean velocity of a particle moving in a random potential in the presence of 
a constant external force. 
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The motion of an overdamped classical particle in a random potential provides an effective description for a variety 
of phenomena, such as the dynamics of dislocations in solids and of domain walls in random field magnets (see Ref. 
|l| and references therein), the relaxation in glasses (see Ref. || and references therein), and the electrical transport in 
disordered solids (see Ref. || and references therein). Recently, this problem has been considered as a phenomenological 
model in the context of glassy dynamics of elastic manifolds in a quenched random mediumou. In this paper we draw 
attention to the mean first-passage time formalism which appears to be very effective for the calculation of different 
characteristics of the random motion. 

We start with the problem of particle diffusion in a long-range correlated random potential and generalize the 
results for the mean squared diffusion amplitude (x 2 (t)} = Aln (t/t r ) obtained in Refs. [?]-§] to the case of nongaussian 
disorder, including the calculation of the exponent (3 and estimates for the prefactor A and the relaxation time t r . 
Second, we investigate the motion of a particle subject to a random potential and driven by a constant external force. 
First, we shall briefly describe the mean first-passage time formalism. 

Consider an overdamped particle in d dimensions subject to the potential V(x) and a gaussian random force rj(t) 
with the correlator (T} a (t)rjp{t')) = 2T5 a p5(t — <'), T denotes the temperature. The boundary conditions involve 
reflecting and absorbing walls S r and S a , respectively. The probability function P(x, t) is obtained as the solution of 
the Fokker-Planck equation 

dP „d 2 P d fdV 



with the boundary conditions P(x € S a ) = and (n • V) P(x E S r ) = 0, with n the unit vector perpendicular to 
the surface S r . The diffusion constant D is related to the amplitude of the random force in the usual way, D = T 
(we set the particle mobility equal to unity). Given the initial condition P(x,i = Q)_=p_^(x — y), the n-th moment 
tn(y) = (*™(y))th °f ^ ne mean absorption time i(y) satisfies the Pontryagin equatior£2nl3 



, d 2 t n _ dVdt 
Qy 2 dy dy 



D ^-iL^ = -ntn-i{y)- (2) 



The above equation has to be solved in a closed region with the boundary conditions t n \g a =0 and (n • V) t n \s r = 0, 
see Refs. [l^ |l2|. Eq. (||) is equivalent to a chain of equations; taking into account that to = 1 we can determine t± and 
proceeding by iteration we find all the moments of t(y) and, consequently, can reconstruct the probability distribution 
function for t(y). Let us apply this formalism to the problem of the ID motion of a particle in a random environment. 

First-passage time moments and distribution function. The equation of motion of a ID overdamped classical particle 
moving in a random potential U(x) takes the form 

dU , , 

, = --+,(*). (3) 
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The Fokker- Planck equation associated with the stochastic equation (|J) is identical to Eq. ([l]) , with D — T, x € 1Z 1 , and 
V(x) = U(x). Let us consider the case where U(x) is a gaussian random potential with correlator ((U(x) — U(y)) 2 ) dis = 
K{x — y), K(u) — > C\u\ , u — > oo. The exponent a is assumed to be positive (if a < 0, the disorder merely leads to 
a renormalization of the diffusion coefficienttHI) . Assume that initially the particle is situated at the point x — y of 
the interval [0, L], with the boundaries and L reflecting and L is absorbing the particle, respectively. Our aim is to 
calculate the disorder-_ayexaged mean first-passage time. The discrete version of the above problem has been studied 
in a number of papersE3ilJ. for the case of a random force with a = 1. 

For the case d = 1 Eq. (||) can be solved exactly, yielding (we remind the reader that y denotes the starting point 
of the particle's diffusion trajectory) 

L yi 

UV) = f[ d yi e u ^' T Je- u toV T t n -i(x Q )dz . (4) 
y o 

Expressing the solution for t n -i(y) in terms of t n -2(y), substituting into Eq. (|j), and proceeding iteratively, we obtain 
the result 

nl f L f y " f L f Vl [ 1 / " \ 1 

tn(v) = J d Vn J dx„_i...y dy x J dxQexpl-\S2u(y l )-^2u{x l )\y (5) 

After averaging over the gaussian disorder we arrive at the expression for the moment (£n(0)) dis 



n 
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(*n( ))dis = 7fk I / dx n _i... / dy x I dx Q 

1 JO JO Jxi JO 

!J2 [K(xi -yj)- K(xi -Xj)- K(y t - yj)] 
2^ y 

where Eq. (^) imposes the In restrictions 

< Xi-x <y t , i = l...n, (7) 

x i < V% < L, i = l...n. (8) 

In general, the integral in Eq. (|J) cannot be calculated exactly. However, in the large distance L limit we can determine 
the integral by the method of steepest descents, thus describing the large-i diffusion. It can be easily seen that the 
integrand reaches its maximum at the point (xi,yi) — (0, L) (note that the restrictions (Q) and (||) are satisfied). As 
K(0) = we obtain 

(*n(0)) dis ~exp(^P). (9) 

The prefactor is determined by the functional dependence close to the saddle-point. Expanding the expression in the 
exponent of (^) and taking into account that K'(0) — we arrive at the final result for the n-th momentli3 of t 

n! / 2T 2 \ 2 " fn 2 K(L) 



The fact that the maximum of the integral in Eq. (||) is realized at the boundary manifests itself through the appearence 
of the first rather than second derivative of K(x) in Eq. (|Io|). 

Using Eq. (0) we can reconstruct the tails of the probability distribution function for the first-passage time t = t((X)- l 
We look for a function of the form P(t) ~ exp (—A In 7 (t/to))- Calculating the moments of P(t) by steepest descentsl 
and comparing with Eq. ( |l0[ ) we find 



'i 
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where io(T,L) is a microscopic time scale (while to accounts for the (dimensional) prefactor in Eq. jlO|), its full 
dependence on L and T cannot be reconstructed from the large-i asymptotics alone). The asymptotic expression (|ll| ) 
is applicable for t > t and produces the strong intermittency observed in the moments (t n ) dis (see Eq. (|U|)). Note 
that the conjecture (in(0)) dis ~ ((ii(0)) dis )™ made in Ref. [TJ| is inconsistent with our findings. 

Large-time diffusion. Using the result (O) we can extract a lot of information concerning the large-t behavior of 
the particle. Suppose that at t = the particle is at the point x = 0. Let us estimate its squared average displacement 
after a time t. The characteristic value of x{t) can be found from the implicit equation (T 2 /K (x)) In 2 (t/to) ~ 1 (see 
Eq. JlH)) defining the characteristic value of x where the probability distribution function P(t) becomes negligible. 
With K(x) ~ C|x| Q , we easily find that 



n2\ 2 / Q 



where t r is a macroscopic diffusion or relaxation time. An estimate for t r is obtained by comparing the characteristic 
barrier Ul ~ V C'L a on scale L with the temperature T. This defines the microscopic diffusion scale 

L T ~(T 2 /C) 1/a (13) 



and its associated diffusion time 



L|/D~ (l/T)(T 2 /C) 2/a . (14) 



Strictly speaking, the problem of the diffusion on a semi-axis which we consider here (the boundary x = is 
reflecting) differs from that of the diffusion on the whole axis. However, the boundary condition at x = affects the 
answer (see Eq. (|l2|)) only by a factor of order unity. 

Eq. (O) can be interpreted in terms of the decay of metastable trapped states. Suppose that a particle leaves 
a metastable state via thermal activation over a barrier of random height. Let the barrier distribution function be 
gaussianCJ, P(U) — (l/\/27rA) exp (— U 2 /2A 2 ). Then the probability distribution of lifetimes i exp(U/T) is given 
by the expression (we assume that t > t , see also Ref. ^0|) 

w-T&A^'k-^-ttM^'M*}- (15) 

o 

Eq. (^||) exhibits the same large-t dependence as Eq. (^lj), implying that the diffusion on the scale L is dominated by 
one deep potential well of characteristic depth A" 1 / 2 (i). This interesting feature of ID diffusion allows us to generalize 
Eq. ( |l2| ) for the case of non-gaussian disorder. 

Assume that the probability for the potential difference U(x + L) — U (x) to be equal to E is given by the function 

P{E) = 2WKV'{Lni/S) CXP { " 2K{L) ) ' (16) 

with K(L) — » CL a , L — > oo. Calculating P(t) in the same way as in Eq. ( |l5| ) above, we arrive at the result (see also 
Ref. H) 

f T 5 \ln(t/t )\ 5 ) 

P W- ex P Jr. >, (17) 



2K(L) 



implying a log i-diffusion of the form 



where t r ~ {^/T) (T s /C) 2 ^ a is the relaxation time for non-gaussian disorder. One can easily verify that Eq. jT^ ) 
implies that In (t n ) dis ~ n' 5 /'' 5-1 '. For 8 — > oo, the barriers in the system vanish and the phenomenon of intermittency 
disappears. 

Eq. ( |l8| ) is consistent with the exactly solvable gaussian random force problem, for which the amplitude before the 
logarithm A = (61/45) T 4 /C 2 is known exactly and the relaxation time t r ~ T 3 /C 2 , see Refs. 0M||. For the case 
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of gaussian disorder (5 — 2) the exponent 4/ a of the logarithm in Eq. ([18]) has been obtained using RG-techiquesa. 
The mean first-passage time method allows us to estimate the amplitudes and characteristic relaxation times and to 



generalize the results to the case of non-gaussian (disorder. Furthermore, Eq. ( JlT[ ) generalizes the results of previous 
investigations of the gaussian random force modelEZI to the case of arbitrary disorder. 

Creep under the action of an external force. Let us apply the mean first-passage time formalism to the problem of 
creep in ID. In a recent papeiu, Le Doussal and Vinokur reported results on the mobility of an overdamped classical 
particle moving in a ID random potential U(x) in the presence of a constant external force / (see also Ref. ||). The 
mean velocity V has been calculated as a function of / and the correlator K(x) of the random potential. The above 
model has been considered as a phenomenological model of glassy dynamics: It turns out that long-range correlations 
of the random potential lead to the glassy response V ~ exp (— 1/T/ M ) as / — > 0, with T the temperature and [i > 
a constant, whereas in the case af short-range correlations F ~ / as / ^ 0. |— ■ 

The problem has been solveda using a generalization of the method introduced by DerridaE^I for discrete models, 
assuming that the random potential is a periodic function of the coordinate x, U(x) = U(x + L). Next, the stationary 
solution P(x) of the Fokker-Planck equation has been found for a fixed current J. Using the conditions P(0) = P(L), 
U(0) = U(L), V = JL, the normalization condition on P{x) 7 and taking the limit L — ► oo, the authors of Ref. |^ arrive 
at the result (see also Ref. 0) 



1 _ 1 

V ~ T 



dse-f s ' T {e^ +s ^- u ^y T ) x , (19) 

o 

where 

L 

(A) x = (1/L) lim / dxA(x). (20) 



For the case of a gaussian random potential with ((U(x) — U(y)) 2 ) dis = K(x — y) the averaging procedure can be 
easily performed, yielding the final result 



1 If, (fx K(x) . 

- = - dxexp(-L r + -±+). (21) 



V T J r \ T 2T 

o 

Let us show how Eqs. (jl^) and ( |2l| ) can be obtained in a simple and elegant way using the first-passage time method. 
This technique has several advantages: i) it does not rely on the periodic continuation of the random potential and, 
consequently, one does not have to worry about the commutation of the two limits t — ► oo and L — > oo; ii) instead 
of solving the stationaxiz-Fokker-Planck equation with a fixed current one can use the well-known solutions of the ID 
Pontryagin equationcjT3 and simply average them over the disorder; Hi) using the first-passage time technique one 
can investigate the finite size effects, the moments of the average velocity distributon function etc., i.e., the information 
obtained is much richer. 

Returning to the ID problem of slu. cwerdamped particle subject to the potential V(x) = U (x) — fx we can solve 
Eq. (S) for the first moment exact ly££rt3, 

L z 

i l(y ) = iy dze v(')/T J dxe~ v ^' T . (22) 



Here we assume that the point L is absorbing, the point a is reflectingcj, and a < y < L. In the limit a — > — oo the 
particle does not feel the left boundary as the external force / is chosen to be positive. Averaging Eq. (|l|) over the 
disorder we obtain 

L z 

(h(y)) dis = ±f dze-W J teef'/T^W-uwv^ (23) 

V -oo 

If the random potential distribution is spatially homogeneous, the average {e^ u ^~ u< ^ x "/ T ) dis is a function of \z — x\ 
and is the same as the translation average we have introduced the new 

variable s = x — z. Eq. (|2^) then takes the form 
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(ti(y)) 



L 



dis 



T 



(24) 



If the integral in Eq. ( p4| ) converges, the mean first passage time averaged over disorder is an extensive quantity. In 
the limit L — y — > oo we can apply the central limit theorem, implying that l/(ti(y)) dis = (l/fi(y)) dis . The mean 
velocity V is simply ((L — y)/ti(y)) dia and we can easily see that Eq. ( pi[ ) becomes equivalent to Eq. (p^). For a 
gaussian random potential, 



(*i(»)>dfa = 



y 



T 



ds exp — 



fs K(s) 
T 2T 2 



(25) 



which is equivalent to Eq. (pl|). 

Relation to discrete models. Let us compare the mean first passage-time averaged over disorder obtained using 
Eq. (22) with that found for a discrete ID random walk in a random force- field. Consider a ID disordered lattice. Let 
us denote by p n the probability of hopping forward n 



The index n enumerates the sites of the lattice. For the random-force problema, 



1 and by q n = 1 ~ p n that of hopping back n 



1. 



exp 


'dF n + 1 ' 
IT 




exp [-^] 


+ exp 


|~dF n + 1 ] 
IT 



(26) 



with d the lattice spacing (note that in this model both time and space are discrete). The force F is a gaussian random 
variable satisfying the conditions (F n F m ) = /j<5„ m /e? and (F n ) = 0. Note that (In (p n /gn)) = 0, see Refs. ^,0. The 
point ft = is supposed to be reflecting and N — L/d is absorbing. In Ref. [15] the estimate (£i(0)) di ~ exp \J3L], 
where /3 = In (q n /p n } = ^i/4T 2 and L — > oo, has been obtained. 

The continuous version of this problem is described by Eq. (||) with a correlator for the random potential given by 
K(x — y) = [i\x — y\. Averaging Eq. (p2) over disorder we obtain 



4T 3 

(*i(y)>dis = ~r 

[I 



exp 



2T 2 y 



eXp l2T^ 



2T 



For y close to the reflecting boundary and for large L the result takes the simple form 

<' (1) (2/)>d, 



4T 3 / uL 



(27) 



(28) 



We find that for both (discrete in space and time and continuous) cases In (ti(0)) di ~ L. The coefficient of propor- 
tionality, however, is different, (3 — /i/2T 2 for the continous model versus (3 = ^/4T 2 for the discrete case. 

We thus have arrived at a different asymptotic behavior for the average first-passage time. The origin of this 
difference is found in the inequivalence of the discrete in space and time and the discrete in space-continuous time 
models: For the latter, consider the lattice master equation describing the probability P n {t) for the particle to be on 
the site n: 



dPn - w P 

—7— — YVn,n+l "n+1 

dt 



W n , n - X P n -i - (W n 



+1,11 



W n -l,n) Pn, 



(29) 



where the hopping probabilities W n;n ±\ and W n ±x,n are determined by the potential!! (as the lattice spacing tends to 
zero, one easily recovers the continuous Fokker-Planck equation (|lf) and thus this model is equivalent to the continuous 
in space and time model). The terms W^n+i-Pn+i and Wra.n-i-Pn-i in Eq. (29) describe the hopping n + 1 — > n and 
n — 1 — ► n, respectively. The contribution (W n ^-\ n + W n —i n ) P n describes the probability to stay at the same site 
n. In a next step let us also discretize the time variable t. Suppose that at time t, P n — 5 nm . In the limit of a finite 
but small time step At the probabilities for the processes, m — > m ± I behave as a\ At and ct2At [a\ and «2 are two 
constants), hence the probability for a particle to stay at the site m after t — > t + At is I — [ax + a^jAt, i.e. the particle 
most likely stays at the same site. On the other hand, for the discrete lattice walks considered in Refs. * nc 

probability of staying at the same site vanishes by definition: The particle may only jump to the neighbouring sites. 
Thus the continuous random walk cannot be obtained as a limiting case of a discrete random walk when both the time 
and the lattice are discrete. 
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Briefly summarizing, we have investigated the motion of an overdamped classical particle in a random potential. 
The moments of the first- passage time averaged over gaussian disorder have been found and are given by Eq. (|ic|). 
The asymptotic form of the first-passage time distribution function has been reconstructed, see Eq. (11), allowing 
us to find the large-i dependence of (x 2 (t)), see Eq. (|T^). The results obtained can be easily interpreted in terms 
of the decay of metastable states. The latter feature has allowed us to generalize the results for (x 2 (t)} to the case 
of non-gaussian disorder, see Eq. (|l^). In addition, we have presented a simple derivation of the mean velocity of a 
particle driven by a constant external force across a disordered medium and have discussed the relation to the discrete 
time random walk on a lattice. 

We thank Stefan Scheidl and Valerii Vinokur for discussions. 
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